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Abstract 

We propose a retrofitted gravity mediation model which alleviates the gravitino 
overproduction from decays of an inflaton and a supersymmetry breaking field. In 
the model, we introduce an approximate £7(1) symmetry under which the supersym- 
metry breaking field is charged, although it is broken by a mass term of messenger 
fields to generate gaugino masses of order the weak scale. In a low-scale inflation 
model, we find regions in which the gravitino overproduction problem is avoided. 



1 Introduction 



Supersymmetry (SUSY) is one of the most plausible candidates for a theory beyond the 
standard model. Since we have not observed any supersymmetric partners of the standard- 
model particles yet, SUSY must be broken in the vacuum. The central issue is how to 
mediate the SUSY breaking effect to the visible sector. Among many scenarios proposed 
so far, gravity mediated SUSY breaking models have been thoroughly and continuously 
studied pp. The gravity mediation has both good and bad points. It can naturally 
generate the \x term of the desired magnitude [21 [3] and may also explain the dark matter 
(DM) abundance by the lightest supersymmetric particle (LSP), while it is plagued with 
the SUSY flavor and CP problems. In spite of its potential problems, the gravity mediation 
has attracted considerable attention as the simplest mediation mechanism of the SUSY 
breaking. 

In recent articles [I] it has been pointed out that there is a new gravitino overproduc- 
tion problem in supergravity (SUGRA). That is, many gravitinos are produced directly 
in the decay of the inflaton if the inflaton has a non-vanishing vacuum-expectation value 
(VEV). The detailed analyses have shown that most of the parameter space in the gravity- 
mediation model of SUSY breaking is excluded by this direct gravitino-pro duct ion process 
together with the gravitino production by particle scatterings in thermal bath. A crucial 
point here is that the gravity mediation has a singlet field Z responsible for the SUSY 
breaking, which mixes with the inflaton field [5j. Because of this mixing the inflaton 
decays into a pair of gravitinos. 

Furthermore, the gravity-mediation model suffers from the Polonyi problem [6]. Since 
the Z field should be completely neutral under any symmetries to generate the gaugino 
masses pQ, the origin of the Z field has no special meaning and hence the value of Z 
during inflation is generically different from the minimum in the true vacuum. Therefore, 
the potential energy of the Z will dominate the universe after the inflation and its decay 
destroys the light nuclei produced by big bang nucleosynthesis (BBN), or if the mass of 
Z is larger than 2 x m 3 / 2 it decays into a pair of gravitinos. (Here, m 3 / 2 is the gravitino 
mass.) The succeeding decay of the gravitino destroys again the light nuclei and ruins the 
success of BBN. It has been recently stressed that this problem is not solved even in the 
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case of dynamical SUSY breaking [7j. 

All above problems are originated from an assumption that the Z field responsible for 
the SUSY breaking is completely neutral. To avoid the problems while keeping the merits 
of the gravity mediation stated at the beginning of this section 0, we propose a gravity- 
mediation model with a Z field charged under some symmetry. In this model the gaugino 
masses vanish at the tree level, since couplings of Z to the kinetic functions of the gauge 
multiplets are forbidden by the symmetry. So we introduce a pair of messengers whose 
mass term breaks the symmetry, to generate the gaugino masses. The one-loop diagrams 
of the messengers, in fact, induce the gaugino masses picking up the symmetry-breaking 
mass term of the messengers. We show that the present gravity-mediation model indeed 
relaxes the gravitino-overproduction problem mentioned above if the inflation scale H m { 
is sufficiently low as H in{ ^ a few x 10 6 GeV. 

The paper is organized as follows. In Sec. [2] we describe the retrofitted gravity medi- 
ation model. We discuss the cosmology of our model in Sec. [31 particularly focusing on 
the gravitino production from both the SUSY breaking field and the inflaton. The last 
section is devoted to conclusions. 

2 A retrofitted gravity-mediation model 

The model is based on a dynamical SUSY-breaking model proposed in [8], which assumes 
an SP(1) gauge theory with 4 chiral superfields Q i (i — 1 — 4) in the SP(1) doublet 
representation, where the gauge index is omitted. Without a superpotential this theory 
possesses a flavor SU(4)p symmetry. We assume, for simplicity, that the flavor symmetry 
is explicitly broken down to an SP(2) F by a superpotential. Thus we introduce 5 gauge 
singlet superfields Z a (a = 1 — 5) and assume the tree-level superpotential 

Wo = X'Z a (QQ) a , (1) 

where {QQ) a denotes a flavor 5-plet of the SP(2) F given by a suitable combination of 
SP(1) gauge invariants Q l Q J . Together with the effective superpotential induced by the 

1 The purpose of this paper is not to solve the SUSY flavor and CP problems, but to alleviate the 
cosmological problems that the original gravity mediation has. 
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strong SP(1) gauge interactions, 

W dyn = X(Pf(QV)-A 4 ), (2) 

the superpotential Wdyn implies that the SP(2) F singlet (QQ) = \(Q l Q z + Q 2 Q 4 ) con- 
densates and we find 

((QQ)) = A 2 . (3) 

We further introduce an SP(2) F singlet superfield Z and consider a tree level super- 
potential, 

W = W + \Z(QQ). (4) 
For the coupling A < 0(1), we find the vacuum, 

((QQ)) = A 2 , ((QQD-O. (5) 

After integrating the massive modes we have the low-energy effective superpotential 

W eS ~ XA 2 Z, (6) 

which yields a dynamical SUSY breaking [8], 

Fz ^ AA 2 . (7) 

Notice that the tree-level superpotential Eq. (jlj) possesses a global U(l) symmetry at 
the classical level, under which the Z fields and Q l transform as 

Z -> e~ l5 Z, Q l -> e + ? s Q\ (8) 

We use this global U(l) to avoid the gravitino-overproduction and the Polonyi problems 
as explained in the introduction q although it is broken by SP(1) instanton effects at the 
quantum level (see also Eq. ([2D). 

In SUGRA the gravitino acquires a SUSY-breaking mass from Eq. (fTI) as pTj 

Fz 

%/2 ^ -y- (9) 



2 A discrete is sufficient for our purpose where the Z, QQ and 'J'S have odd parity of the Zi 
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QQ 
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U(1)r 





+2 


+2 





U(l) 


+1 


-1 


-1 


+1 



Table 1: The charges of U(1)r and U(l). 



Here and in what follows, we adopt the Planck unit, Mp = 1 unless otherwise stated, 
where Mp ~ 2.4 x 10 18 GeV is the reduced Planck scale. For a generic Kahler potential 
squarks, sleptons and Higgs bosons acquire the SUSY-breaking soft masses of 0(m 3 / 2 ). 
However, the gauginos in the SUSY standard model (SSM) remain massless [HI [10] , since 
the Z does not have couplings to the gauge kinetic functions. In fact, the interaction, 

J d 2 9ZWX, (10) 

is forbidden by the U(l) symmetry in Eq. (jHJ), where W£ are chiral superfields for gauge 
multiplets. Therefore, we need a breaking term of the global U(l) symmetry to generate 
gaugino masses. Otherwise, the dominant contribution to the gaugino masses comes only 
from the scale-invariance anomalies at the quantum level [TT], which may be of order 
x 7713/2- For m 3/ / 2 = 100 GeV — 10 TeV we have the gaugino masses of order 1 — 100 
GeV which is excluded already by experiments. Thus, to have larger gaugino masses 
we introduce a pair of messengers \I/ and ^ whose mass term breaks the global £7(1), 
assuming that they transform as 5 and 5* of SU (5) gut, respectively. 
Then, the Z field has a superpotential with the messenger fields as 

^ mesS enger = kZV* + MM. (11) 

The U(l) and U(1)r charges for relevant superfields are given in Table 1. We see that 
the messenger mass term M breaks the global U(l) symmetry. Here and in what follows, 
we assume that, at the breaking scale M, the effect of the U(l) breaking appears only in 
the mass term of \I/ and ^ Q The integration of the messengers give rise to the gaugino 

3 If one allows any [/(l)-breaking operators suppressed by powers of M, there are such dangerous 
operators as K ~ M* Z\<f>\ 2 +h.c. which induce the severe gravitino overproduction, where <f> denotes the 
inflaton field. 
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masses as p2] @ 

for < = w- ( 12 ) 

Here, 777,1,2,3 and 0:1,2,3 are the gaugino masses and the gauge coupling constants for 
U(l), SU(2) and ^[/(S) in the SSM. We have used the SU(5)gut normalization for the 

£7(1) gauge coupling constant. For m 3 ~ 1 TeV we have 

h W 

— ~ 10 5 GeV. (13) 

Notice that the global U(l) and U(1)r charges for the operator QQ are the same as 
the tyty and hence the dynamical quarks may naturally have a mass term M'QQ with 
M' ~ O(M). In the text, we have redefined the Z field by a shift , Z — > Z — M', to absorb 
the mass term for the dynamical quarks Q. However, this shift of the field induces, for 
instance, a linear term of Z in the Kahler potential (see Eq. (12"T|) ). In the following, we 
adopt the origin of Z as that obtained after the shift. 

We should mention here that the Giudice-Masiero(GM) mechanism [2] for generating 
a SUSY- invariant mass term (called the /x term), fiHH, for Higgs doublets does not work, 
provided that the Higgs multiplets, H and H are neutral of the global U(l) symmetry ^. 
However, they receive the /i term though the following superpotential [3], 

W = C(l + hHH), (14) 

where C = is the constant introduced to cancel the vacuum-energy density \Fz\ 2 for 
the SUSY-breaking. Here, we have assumed that HH carries a vanishing U{1)r charge. 
Then, we find the // parameter for the Higgs mass as 

A* = ms/2 x h. (15) 

In the following discussion we restrict ourselves to the parameter region of so-called 
gravity mediation, that is, m 3 / 2 ~ lOOGeV — 10 TeV and the gluino mass m 3 ~ 1 TeV. 
This implies from Eqs. (|9D. (fT2|) and ([TBI 

JFz = ^AA ~2x 10 10 GeV - 2 x 10 n GeV, (16) 



4 Although the scalar trilinear couplings are suppressed at the dynamical scale in this model, it is 
possible to induce sizable contributions by introducing the Yukawa interactions, YsSHty + YgSH i &, with 
the SM singlet S, assigning the U(l) and U(1)r charges for S and H, H properly [T31 W£\ . 

5 If one assumes the [/(l)-charge +1 for HH, the Kahler coupling HH is allowed and the GM 
mechanism works. 
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and 

A ~ 3 x lO-^-^GeV" 1 . (17) 
M y ' 

We see that the SUSY-preserving vacuum, (Z) = -M/k ~4x 10 15-17 GeV with (\P 1 3>) ~ 
—XA 2 /k, is far from the dynamical scale \Z\ ~ A in Eq. ffl6|) . and hence the SUSY- 
breaking vacuum, \{Z) \ ^ A and Fz — AA 2 , is practically stable. Indeed, the messenger 
fields are not tachyonic at the origin Z = and the tunneling rate is also suppressed as 
long as M> ykXA, which is satisfied unless k is extremely small |15j . 

3 Cosmology 
3.1 Polony i problem 

Let us first discuss the Polonyi problem in the present model. We assume that the SP(1) 
hadrons have masses of order An A and hence above the scale \Z\ > \Z*\ = AnA/X, the 
Polonyi field Z does not receive effects from the SP(1) strong interactions. The potential 
of Z is therefore very flat above \Z*\. On the other hand, for \Z\ < \Z*\ the Z acquires a 
larger SUSY-breaking soft mass from the loop diagrams of the SP(1) hadrons, and it is 
given by 

m z ~ -^A 3 A, (18) 

l07T^ 

where r] is a numerical coefficient which is expected to be order unity. The potential for 
Z can be approximated by [IE] 



£ rn\ j2 for \Z\ > \Z* 



m 2 z \Z\ 2 for \Z\ < \Z, 
Yl(Z)~{ 2 . (19) 

l 3/2 



where we have set the cosmological constant to zero at the origin, and £ ^ 0(1) is a real 
constant. Note that, for \Z\ > \Z*\, the potential is nonzero due to the SP(1) gaugino 
condensation, and the curvature of the potential is given by |7j 

Vl'(Z),-^ + OK /2 ), (20) 

where the first term comes from the perturbative wavefunction renormalization of Z, while 
the second term represents the contribution from the gravity mediation. For A ~ 0(1) 
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Figure 1: Schematic potential of the scalar component of the Z field. 



and \Z\ <C 1, the first term dominates over the second term, and so, we will focus on the 
first term from here on. A schematic potential for the Z boson is shown in Fig. [U 

If the global U(l) symmetry discussed in the previous section is exact, the Kahler 
potential for the Z is a function of Z^Z and the origin of the Z is most likely the poten- 
tial minimum not only at the present, but also during the inflation. However, since we 
introduce a mass parameter M for the symmetry breaking, it is natural to consider the 
Kahler potential contains breaking terms. Thus we consider the Kahler potential for the 
Z as@ 

K= (5M ] Z + h.c.) + ZZ ] + ■■■, (21) 

where 5 is a constant of 0(1). 

We see that the minimum of the Z potential shifts from that at the vacuum during 
the inflation. The effective potential for Z during inflation is given by ^ 

V(Z) ~ e K (3Hl { )+V L (Z), 

~ 3#f nf (\Z\ 2 + 5M^Z + 5MZ^ + ■■■) +V L (Z). (22) 

If the Hubble parameter during the inflation satisfies H m f > mz, the shift is therefore [j] 

\AZ\ ~ \8M\ > Z*. (23) 



6 Even if the Kahler potential does not contain such symmetry breaking term from the beginning, it 
is radiatively induced as K ~ —L^kM^Z + h.c. 

7 It should be noted that, during inflation, Vl{Z) is absent for a high-scale inflation model with 
H m i > 4nA. However, this does not change our argument. 

8 Such a large deviation from the origin is dangerous not only because it induces the cosmological 
disaster discussed below, but also because the Z field may roll down to the supersymmetric vacuum after 
inflation. 
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After inflation, the Z field starts to oscillate with an initial amplitude |AZ|, when the 
Hubble parameter becomes comparable to the curvature of the potential at Z ~ \AZ\. 
Since the potential has an approximate £7(1) symmetry and is almost flat for \Z\ > \ZJ\, 
the Z field experiences spatial instabilities and soon deforms into Q-balls [XT] with a 
typical charge [18], 

IAZI 4 



-4 



Q = 6x10 \v£(Az)r 

where we have used f|T2|) in the second equality. After the Q-balls are formed, the energy 
density of Z decreases as a -3 (a is the scale factor) like a non-relativistic matter. Since 
the lifetime of the Q-ball is rather long as shown below, the Z field will dominate the 
energy of the universe before the decay. 

How does the Z field in the form of the Q-balls decay? The Z field decays into a pair 
of the gravitinos as well as the SM fields. Since the Z decays only through (effective) 
higher dimensional operators, the decay rate of the Q-balls is unlikely to saturate the 
geometrical upper bound [IS]. For the moment we assume this is not the case. Then Z 
decays in the entire volume of the Q-balls. Inside the Q-balls, the Z field rotates in its 
internal space with a constant angular velocity, and the Q-ball solution is of the form: 

Z(r,t) = Z(r)e iujt . (25) 

Here uj ~ ^"(AZ)! 1 / 2 is the angular velocity and Z(r) represents the radial profile of the 
Q-ball with a radial coordinate r, ranging from r = to r = Rq, where Rq is the radius 
of the Q-ball. The field value at the center r = 0, Z(0), is roughly equal to AZ. For later 
use, let us express to in terms of the gluino mass as 

^ £ {m) m: " (26) 

where we have used ( 1121) and (12"U|) . The classical-field configuration (|25|) is interpreted as a 
condensate of the Z-particles with energy uj per unit quanta and with a macroscopic num- 
ber density ~ u\AZ \ Therefore one can use the perturbative decay rate by substituting 
uj for the mass of Z 



9 Since the finite size of the Q-ball implies a finite momenta, p^SRq 1 ~ u), for a quanta in the 
condensate, numerical factors may change by O(l) in the following decay rates. 
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Since the Z field has a nonzero F-term \Fz\ ~ ^3/2 even for |Z| > \Z*\ due to the 
gaugino condensation, it decays into a pair of the gravitinos. The decay rate is given by 

1 u 5 



T z (Z^2^ 3/2 ) ___ , . 

' 967T 777.3,2 



where m 3 and 7773/2 are evaluated at the weak scale. In the goldstino picture, the above 
decay into the gravitinos is induced by the four point interaction \Z\ A in the Kahler 
potential, which arises from the wavefunction renormalization of Z. 

In addition, Z interacts with the SM gauge sector via the messenger loops. Since 
uj <C M, we can write an effective coupling as 

£ = -/^^« , + h - c '> < 28 > 



where Z is the superfield. This interaction is expanded as 

7„ ,,1 

+ h.c, (29) 



£ — — a — vjZ 

4ttM 



where we neglected terms with higher orders of k(Z)/M. Then the decay rate into the 
SU (3)c gluons is 

k 2 ( a-A 2 uj 3 



~ ntt/ M 2 ' 



.(Jib 

while that into the gluinos is 



* 6xlO-GeV^j^ 

where we have assumed uj is much larger than 27773. 

The Z also decays into the SM matters via the top Yukawa coupling, since it has a 
linear term in the Kahler potential [20]. However, the rate is smaller than that into the 
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gauge bosons. The decay via anomalies in SUGRA [23] is also suppressed by a loop factor, 
so we neglect them here. 

One can check that the decay rates (T27j) and (|30|) (or (I3T]) ) are much smaller than the 
geometrical upper bound on the Q-ball decay rate: [19] 



1 



dQ 



dt 



The upper bound can be thought of as a dissipation rate of relativistic decay products. If 
the decay products are fermions (e.g. the gravitinos and gauginos in our case), and if the 
perturbative decay rates obtained in a way described above exceed the bound, the decay 
processes inside the Q-ball would be suppressed by the Pauli blocking since the decayed 
products would fill the phase space. Then, the dissipation rate would determine the decay 
rate of the Q-ball. In our case, however, since the decay proceeds only through the higher 
dimensional operators and the perturbative decay rates are so small, such suppression 
is absent. On the other hand, if the decay products are bosons (e.g. the gluons in our 
case), there is no Pauli blocking inside the Q-ball. However, if the bosons acquire a large 
mass due to interactions with the scalar field that form the Q-ball, the decay inside the 
Q-ball might be kinematically blocked. Then the decay rate of the Q-ball would become 
again Fq. In the present case, however, since the Z field is singlet under the SM gauge 
symmetry, the gluons are massless inside the Q-balls. Thus the decay rate of the Q-ball 
is given by the perturbative decay rate. 

It is illustrative to take the ratios of the above decay rates: 

\\ (33) 



T Z (Z -> 2^3/2) _L(_^ 



Therefore, for (X/kS) 3> 1, the dominant decay mode is that into the gravitinos, while the 
gluino production dominates over the others for 0.1 <{X/k5) < 0(1). For (X/k6) < 0.1, 
the decay into the gluinos is kinematically forbidden (see ( 1261) ). and the gluon production 
dominates over the other two channels. Note that, for (X/k5) ^> 1, the effective mass of 
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Z, uj, is much larger than the weak scale and the decay into a pair of the gravitinos is 
kinematically allowed. 

First let us consider the case of (X/k5) ^> 1, in which the gravitino production is the 
main decay mode. The total decay rate of the Z field is given by r^ ota1 ^ ~ r z{Z — > 2^ 3 / 2 ). 
Then, after the decay of Z, the universe is dominated by the gravitinos with a small 
amount of the entropy produced by the decay into the gluons and the gluinos. The 
gravitino-to-entorpy ratio is 



i 

r2, 



Y (Z) _ 2g 3/2 3 hi 2 gA 4 r^Ii) 



3 / 2 " (l-B 3/2 ) 3 /M W \10 J yLz ' 

» gxio-^^V'f-V 1 ^ 1 ^)" 1 ^ (35) 

V10.757 V0.17 VTeVy VTeW ' v ; 

where counts the relativistic degrees of freedom, and _B 3 / 2 — 1 denotes the branch- 
ing ratio of the gravitino production. Thus the gravitino abundance is too large to be 
compatible with the constraints from BBN which range from O(10 -16 ) to O(10~ 14 ) for 
m 3/2 = 100 GeV - lOTeV [2H 122] (see also (jSgjl). 

On the other hand, the gluino production dominates for 0.1 ^(X/k8) ^ 0(1). Since the 
decay temperature is rather low, the resultant lightest SUSY particles (LSP) will easily 
overdose the universe. In order to avoid the overproduction of the LSP, the effective mass 
u should be smaller than 2m 3 , i.e., (X/k5) i$ 0.1. Using (I3"0l as the total decay rate, the 
decay temperature becomes 



_2. 



Td s *m 



V10.757 \6k) vo.i; VTeW VTeW ' 

BkBvf-^VV— yf^" 1 - (36) 
V 10.757 VTeW VTeW y J 



< 



Since the gaugino mass is proportional to the gravitino mass, the decay temperature in- 
creases as the gravitino mass becomes larger. However, even for m 3 / 2 — ^3 = O(10)TeV, 
it is still smaller by several orders of magnitude than the lower bound on the decay tem- 
perature by BBN [2TI E2] , the cosmic microwave background and the large scale struc- 
ture [23J. Thus we conclude that once Z field deviates from the origin during inflation 
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by AZ ~ 8M, the Polonyi problem associated with the Z field spoils the success of the 
standard cosmology. 

Now let us turn to a low-scale inflation scenario, satisfying H in { < mz- If the initial 
position of Z is beyond Z*, it will settle down at \Z\ ~ \8M\ during inflation, and the 
Polonyi problem jeopardizes the scenario as before. In the following we assume that this 
is not the case. The shift of the minimum is then given by 

H 2 

\AZ\ ~ -if\8M\. (37) 
We require \AZ\ < \Z*\ to avoid the above Polonyi problem, which leads to a constraint 
H M < 2 x 10 6 GeV^L (H*) 1 (^) * . ( 38 ) 



After inflation, the Z starts to oscillate with an initial amplitude AZ given by (137]) . and 
soon decays into a pair of the gravitinos, since the rate is enhanced especially when mz 
is much larger than 777,3/2 (see the equation above (1271) ). The abundance of the gravitinos 
produced by the Z decay is 



where Tr denotes the reheating temperature of the inflaton, and we have used ( 1381) for the 
last inequality. We find that the gravitino abundance is much smaller than the previous 
case (135]) . and it can be compatible with the BBN bounds. 

To sum up, the Polonyi problem associated with the Z field excludes high-scale infla- 
tion models, and only low-scale inflation models satisfying (138"]) may be able to circumvent 
the problem. 

3.2 Gravitino production from Inflaton decay 

The gravitino production from the inflaton decay is a quite generic phenomenon. In fact, 
for a large-scale inflation model, the inflaton directly decays into the SUSY breaking 
sector [201 El], producing the gravitinos, while the gravitino pair production becomes 
effective for a low-scale inflation model. 
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Let us first consider inflation models with an inflaton mass larger than the dynam- 
ical scale An A, which is typically the case for large-scale inflation like a hybrid inflation 
model [25]. The inflaton then decays into the SUSY breaking sector through the following 
processes. As pointed out in Ref. [20J, the inflaton decays via the Yukawa coupling (@J, 
producing the scalar and fermionic components of Z and the hidden (s) quarks Q % . Note 
that the fermionic components of Z is the goldstino which will be eaten by the gravitino, 
and that the scalar Z dominantly decays into a pair of the gravitinos. In addition, the 
inflaton decays into the SP(1) gauge sector via the anomalies in SUGRA [24]. The pro- 
duced gauge bosons/gauginos form jets, producing the SP(1) hadrons. In the decays of 
the SP(1) hadrons, the gravitinos are produced. The detailed analyses [201 El] actually 
show that the gravitino production through the above processes excludes almost entire 
parameter spaces for the large-scale inflation models. Therefore, even if one circumvents 
the Polonyi problem discussed in the previous section (e.g. by fine-tuning the linear term 
of Z in the Kahler potential), the high-scale inflation models still suffer from the severe 
gravitino overproduction problem. 

In the rest of this section, we focus on low-scale inflation models with the inflaton 
mass smaller than 47rA. Then the spontaneous decays [20, El] do not occur, since the 
SUSY breaking fields typically have a mass of the dynamical scale. Instead, we need 
to take account of the gravitino pair production from the inflaton. The gravitino pair 
production occurs even when the inflaton has the minimal Kahler potential. When the 
Kahler potential is minimal for the inflaton field, the gravitino pair production rate is [5] 



where (0) denotes the VEV of the inflaton, and the upper index (0) is to remind us 
that the rate is for the minimal Kahler potential. Although the gravitino production is 
suppressed if ^> m z , one cannot expect too large hierarchy between m$ and m z , 
because we are considering low-scale inflation models with < 47rA. 

Now we consider non-renormalizable couplings of the inflaton with Z. Actually, those 
couplings that induce the mixings between and Z in SUGRA enhance the gravitino 



rJ )) (0-2^/ 2 ) * ^<0>XxF(°) 



(41) 



with 




(niz/m^Y for m^^mz 
1 for 77i0 <C mz 



(42) 
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production. The relevant mixings arise from the following interactions: 



K = P^Z + ^l^ZZ + h.c. 



(43) 



where flip) are numerical coefficients. The presence of those interactions is rather generic, 
since even though we have assumed that such t/(l)-breaking operators are absent at the 
cutoff scale, these are radiatively induced during the evolution running down to lower 
energy scale through the U(l)-breaking operator in Eq.t jTTT) . At the inflaton mass scale 
-C M, /?i(2) are estimated as 



0i 

02 



^M*ln^-. 

16tt 2 |M| 2 ' 

5C (kM*) 2 



16tt 2 \M\ 2 ' 

where we have introduced a non-renormalizable interaction, 

K = C|0| 2 |^| 2 , 



(44) 
(45) 



(46) 



with a numerical coefficient ( of order unity. The decay rate into a pair of the gravitinos 
is0 



with 



f |/3i| 2 m 4 



F W = 



327T 

z +|/? 2 | 2 for m^>m z 



(47) 



3 m l/2 m l 



W ™J , | fl |2 



(48) 



» + |&| — fo r m^^niz 
3 m^, 2 \m z / 



where we have neglected interference terms for simplicity. 

Combining (j4TT) and ( 1471) . the gravitino abundance becomes 

2^3/2) 3T« 



3/2 



-1 (total) 



v2(W lf0 15 Gevj vl0 10 GeW vl0 6 GeW ' 1 } 
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where F\ is the total decay rate of the inflaton, and is related to the reheating tem- 
perature Tr as 

r (total) _ / 9* \ 2 T 2 /c- n \ 

L * ~[~w) { } 

We have also defined F = T (0) + T (1) by ignoring the interferences. 

Let us compare the result ( )49j) to the gravitino production from the Z field associated 
with the Polonyi problem. We notice that the dependence of on the reheating tem- 
perature is different from that of Kj/2 i n (1391) . Importantly, the gravitino overproduction 
problem cannot be solved simply by reducing the reheating temperature, since Yyl is 
inversely proportional to Tr. This makes it nontrivial whether there exist cosmo logically 
allowed parameter regions. 

3.3 Example 

Let us consider a new inflation model [26] as an example of the low-scale inflation models. 
The Kahler potential and superpotential of the inflaton sector are written as 

K(<f>,^) = |0| 2 + ^|0| 4 , 

WW = v 2 4>-^r +1 - (si) 

where the observed density fluctuations are explained for v ~ 4 x 10~ 7 (0.1/g) 1//2 and 
/t ^$0.03 in the case of n = 4. We assume n = 4 in the following, since the Hubble 
parameter during inflation likely exceeds the bound (|38|) for n > 4 [27J. After inflation, 
the inflaton <p takes the expectation value (0) ~ [y 2 / g) l ^ n . The inflaton mass is given by 

~ nv 2 / ((f)), and the gravitino mass is related to v as m 3 / 2 — nv 2 (0) /{n+ 1), since 
the inflaton VEV induces the spontaneous breaking of the i?-symmetry, namely a nonzero 
(W) . Precisely speaking, v has a weak dependence on Tr via an e- folding number and on 
k, and (0) depends on these parameters as well. Therefore one can express the inflation 
scale v and the coupling g as functions of Tr and k, once the WMAP normalization 
of the density fluctuations |28j is applied. In the following numerical analyses, we take 
into account these corrections. 

In the numerical analyses, we estimate the gravitino abundance. In addition to the 
non-thermally produced gravitinos, (|39j) and (|49|) . we also include the contribution from 
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the thermally produced gravitinos: [29jl£ 
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(52) 



The gravitino abundance is severely constrained by BBN as [31 
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(53) 



for the unstable gravitino with a hadronic branching ratio Bf, — 1. 

We show contours of the gravitino abundance Y 3 / 2 for the new inflation model (thin 
solid (red) lines) in Fig. El for several sets of (A, k) = (0.5, 0.1), (1, 0.1), (0.5, 0.01) and 
(1,0.01) labeled by (A), (B), (C) and (D), respectively. We also show the parameter 
space consistent with the BBN bounds ( I53p . which is enclosed by the thick solid (green) 
lines. We have chosen the other parameters of the SUSY breaking sector as 5 = 0.1 and 
C = 1. We have taken rj — 1, which should be in principle determined by the strong 
dynamics, and m 3 = ITeV as a reference value. The parameters of the new inflation 
model, k and n, are chosen as k = 10 -2 and n = 4, while g is not an independent 
parameter and is determined by the other parameters. The reheating temperature Tr is 
regarded as a free parameter, by assuming proper couplings of the inflaton with the SM 
sector [27], though the spontaneous decays via the SUGRA effects [20] provides the lowest 
reheating temperature as Tr ^ 0(1) GeV. 

From Fig. [21 one can see that the gravitino abundance becomes larger for the heavier 
gravitino. It also tends to increase when Tr is both raised and lowered. These behaviors 
are mainly due to the non-thermal productions of the gravitinos, (1391) and (|49p . and it 
can be understood as follows. For larger 7713/2 and Tr, the gravitino production from the 
SUSY breaking field becomes important, while the gravitino from an inflaton dominates 
for smaller Tr. Here it should be noted that the inflaton mass is positively correlated 
with the gravitino mass in the model (|51|) . 



10 The gluon in the hot plasma might decays into the gravitino due to the thermal corrections |30j . 
This, however, changes the gravitino production rate only by a factor. 
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Since the gravitino abundance is sensitive to the model parameters, the allowed regions 
change for a different set of (A, k). When A increases, the gravitino production from the 
SUSY breaking field (|39|) is suppressed, while more gravitinos are produced by the inflaton 
decay (see (1491) ). Thus the allowed region shifts upwards as can be seen by comparing the 
panel (B) with (A) (or (D) with (C)) in Fig.[2j On the other hand, as k becomes smaller, 
the abundances of the gravitino produced both from the Z and the inflaton decrease (see 
(1391) and (SHj)). Comparing the panel (C) with (A) (or (D) with (B)) in Fig. El we find 
much broader allowed region for smaller k. 

We notice that for smaller m 3 / 2 and relatively large Tr, the contours of the grav- 
itino abundance tends to become independent of m 3 / 2 , especially when the non-thermal 
gravitino production is suppressed. This means that the abundance is dominated by the 
thermal production, and therefore it is determined solely by Tr there. 

One may be interested in the regions of higher reheating temperature such as Tr ^ 10 6 
GeV, where non-thermal leptogenesis [32| 133} 134"] may be able to explain the baryon 
asymmetry of the universe. For the panel (C) in Fig. [2j we find allowed regions with 
T R > 10 6 GeV for m 3/2 > 7 TeV and m 3/2 = 100-400 GeV, and similarly, for m 3/2 > 4 TeV 
and m 3 / 2 = 100 — 500 GeV in the panel (D). Notice that the thermal production of 
the gravitino (|53|) imposes constraints on the reheating temperature as Tr ^ a few x 
10 6 GeV for m 3 / 2 = 100 — 500 GeV, while it becomes significantly relaxed for m 3 / 2 ^ 4 TeV. 
Therefore, in the panel (D), the upper bound on the reheating temperature is almost the 
same as that from the thermal gravitino production. 

We have varied A, k and 5 to see how the constraints depend on the parameters in the 
SUSY breaking sector. In Fig. [31 we have examined whether or not the allowed region 
exits for m 3/2 = 100 GeV - 10 TeV and Tr = 1 GeV - 10 10 GeV (left panel), and for 
m 3/2 = 100 GeV - 10 TeV and Tr = 10 6 GeV - 10 10 GeV (right panel). One can see that 
A is bounded below for fixed k, since too small mz makes the Polonyi problem worse 
(see ( 1391) ). Although not explicitly shown in the figure, A cannot be too large since more 
gravitinos are produced from the inflaton. For smaller k and 5, it becomes easier to satisfy 
the BBN bounds, though too small k upsets the stability of the SUSY breaking vacuum. 

As a comment, although we have taken account of the LSP production from the 
gravitinos, it does not give any meaningful constraints on the parameters in which we are 
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interested lj. Indeed, the abundance of the LSP produced from the gravitino is always 
negligible for the LSP mass ^0(l)TeV, as long as we require the BBN constraints on 
Y3/2 ( |53|) to be satisfied. 

So far, we have considered the new inflation model (ED), in which the inflaton mass 
and the gravitino mass are correlated. The relation does not hold if we consider a two-field 
new inflation model [35]. For instance, even for larger gravitino mass, we can take the 
inflaton mass lower than the model ( 15T1) . Then the gravitino production from both the 
Polonyi field and the inflaton can be suppressed. 



4 Conclusions 

The gravity mediation provides a simple way to mediate the SUSY breaking to the visible 
sector, and so, it has been one of the main target of research. However, recent observations 
on the gravitino overproduction from an inflaton and the Polonyi problem revealed that 
the gravity mediation is faced with the cosmological embarrassment, which drives the 
scenario into a corner, hinting that some improvement is needed. In this paper, we 
have proposed the retrofitted gravity mediation model to circumvent the cosmological 
problems. We have introduced an approximate U(l) symmetry under which the SUSY 
breaking field is charged to alleviate the gravitino overproduction from an inflaton and a 
supersymmetry breaking field. Indeed, we have found such regions for a low-scale inflation 
model that all the superparticles, especially, the gauginos as well as the gravitino, have a 
mass around the weak scale and the cosmological bounds on the gravitino abundance are 
satisfied. 

Specifically, there are allowed regions that the reheating temperature is larger than 
10 6 GeV for m 3 / 2 ^ 4TeV and m 3 / 2 = 100 — 500GeV. These regions are attractive since 
non-thermal leptogenesis may be able to explain the baryon asymmetry of the universe. 
For the gravitino mass heavier than 4TeV, the squark and sleptons acquire large masses 
of 0(m 3 / 2 ), and therefore the problems of large flavor changing neutral currents and CP 
violation become mild. Further, if the gaugino (and/or higgsino) masses are so light as 

11 Here we have neglected the thermal production in the estimation of the LSP abundance. It is highly 
dependent on the mass spectrum of the visible sector. 
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Figure 2: Contours of the gravitino abundance Y3/2 and the BBN constraints on the 
(777.3/2, Tr) plane, in the case of the new inflation model. The regions surrounded by the 
thick solid (green) lines are consistent with the BBN constraints ([53]) . We have set the pa- 
rameters as: 8 = 0.1 and k = 10" 2 , varying (A, k) as (A, k) = (0.5, 0.1), (1, 0.1), (0.5, 0.01) 
and (1,0.01) labeled by (A), (B), (C) and (D), respectively. 



O(100) GeV, the particle spectrum resembles that in the focus-point region [36J, and the 
lightest neutralino may be able to explain the current DM abundance [37] . 

Finally let us comment on the possible extension of the model. From the numerical 
analyses, we have seen that smaller k and 8 are cosmologically favored since the gravitino 
overproduction problem gets greatly relaxed. Such a suppression of k and 8 can be realized 
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Figure 3: BBN constraints on A and k for 5 = 0.01,0.1 and 1 (from left to right). The 
regions below the lines are allowed, in the sense that there exists a set of (777.3/2, Tr) 
consistent with the BBN results in the range, m 3/2 G [10 2 , 10 4 ]GeV and T R e [1, 10 10 ]GeV 
(left panel); m 3/2 E [10 2 , 10 4 ]GeV and T R e [10 6 , 10 10 ]GeV (right panel). 

easily by imposing the Z2 symmetry. In fact, if Z and QQ are odd under Z2, the Yukawa 
coupling k and the mass M' should break the Z 2 symmetry. Thus k and 5 = M'/M 
become naturally small. However it is noticed that k cannot be too small, since the 
stability of the SUSY breaking vacuum would be upset due to the tachyonic messengers 
and/or the decay via the tunneling to the SUSY preserving vacuum. 
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